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CHEMICAL  KINETICS  OF  HYDROCARBON  IGNITION  IN  PRACTICAL 

COMBUSTION  SYSTEMS 
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Lawrence  Livermore  National  Laboratory 
Livermore,  CA  94550,  USA 


Chemical  kinetic  factors  of  hydrocarbon  oxidation  are  examined  in  a  variety  of  ignition  problems.  Ignition 
is  related  to  the  presence  of  a  dominant  chain-branching  reaction  mechanism  that  can  drive  a  chemical 
system  to  completion  in  a  very  short  period  of  time.  Ignition  in  laboratory  environments  is  studied  for 
problems  including  shock  tubes  and  rapid  compression  machines.  Modeling  of  the  laboratory  systems  is 
used  to  develop  kinetic  models  that  can  be  used  to  analyze  ignition  in  practical  systems.  Two  major  chain¬ 
branching  regimes  are  identified,  one  consisting  of  high  temperature  ignition  with  a  chain  branching 
reaction  mechanism  based  on  the  reaction  between  atomic  hydrogen  with  molecular  oxygen,  and  the 
second  based  on  an  intermediate  temperature  thermal  decomposition  of  hydrogen  peroxide.  Kinetic  mod¬ 
els  are  then  used  to  describe  ignition  in  practical  combustion  environments,  including  detonations  and 
pulse  combustors  for  high  temperature  ignition  and  engine  knock  and  diesel  ignition  for  intermediate 
temperature  ignition.  The  final  example  of  ignition  in  a  practical  environment  is  homogeneous  charge, 
compression  ignition  (HCCI),  which  is  shown  to  be  a  problem  dominated  by  the  kinetics  of  intermediate 
temperature  hydrocarbon  ignition.  Model  results  show  why  high  hydrocarbon  and  CO  emissions  are  in¬ 
evitable  in  HCCI  combustion.  The  conclusion  of  this  study  is  that  the  kinetics  of  hydrocarbon  ignition  are 
actually  quite  simple,  since  only  one  or  two  elementary  reactions  are  dominant.  However,  many  combustion 
factors  can  influence  these  two  major  reactions,  and  these  are  the  features  that  vary  from  one  practical 
system  to  another. 


Introduction 

In  many  practical  steady  combustion  systems,  ig¬ 
nition  is  simply  a  means  of  starting  the  system  on  its 
way  to  steady  state.  Performance  and  emissions  are 
essentially  independent  of  ignition  in  such  systems 
as  boilers,  furnaces,  and  burners.  However,  in  other 
practical  problems,  ignition  has  a  great  influence  on 
performance,  emissions,  and  other  characteristics, 
and  ignition  can  explain  the  performance  of  the  en¬ 
tire  system. 

Ignition  can  depend  on  physical,  chemical,  and 
mixing  and  transport  features  of  a  problem  and  in 
some  cases  on  heterogeneous  phenomena.  Excellent 
reviews  of  ignition  can  be  found  in  current  sources 
[1-4]  describing  thermal  feedback,  chemical  kinetic 
chain-branching  reactions,  and  other  elements. 
However,  ignition  in  general  is  an  enormous  subject, 
and  the  present  work  cannot  provide  a  thorough 
treatment. 

This  paper  focuses  on  chemical  kinetic  factors  in 
practical  systems,  with  special  attention  on  ignition 
in  automotive  engines.  Recent  advances  in  experi¬ 
mental  and  kinetic  modeling  capabilities  have  pro¬ 
vided  new  insights  into  ignition,  offering  new  pos¬ 
sibilities  for  control  strategies  and  new  classes  of 
combustors. 


General  Features  of  Ignition 

Many  interesting  systems  involve  chain  reactions, 
such  as  nuclear  reactors,  with  neutrons  as  chain  car¬ 
riers.  Chain  termination  is  provided  by  neutron  ab¬ 
sorbers,  and  chain  branching  is  associated  with  the 
term  supercritical,  normally  a  condition  to  be 
avoided.  Populations  of  living  organisms  obey  the 
same  reactivity  laws.  Organisms  can  grow  exponen¬ 
tially  via  chain  branching  until  limited  availability  of 
nutrients  first  stabilizes  the  population  and  eventu¬ 
ally  quenches  the  system  via  chain  termination. 

Reacting  systems  follow  the  general  equation 


where  n(t)  is  the  number  of  chain  carriers.  This  in¬ 
dicates  that  the  change  in  number  of  neutrons  in  a 
reactor  depends  on  the  number  of  neutrons  available 
to  produce  further  neutrons  or  the  change  in  the 
number  of  bacteria  depends  on  the  number  of  bac¬ 
teria  which  can  reproduce,  a  reflection  that  the  cur¬ 
rent  population  of  chain  carriers  produces  the  next 
generation  of  chain  carriers. 

Solution  of  equation  1  leads  to  the  expression 

n(t)  —  n( 0)  exp(kf)  (2) 

When  k  <  0,  the  result  is  exponential  decay;  when 
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k  >  0,  the  radical  population  experiences  exponen¬ 
tial  growth.  Thus,  k  =  0  is  equivalent  to  criticality  or 
steady  combustion  (or  population  stability).  In  a 
chemically  reactive  system,  the  coefficient  k  is  an 
average  over  all  reactions  taking  place,  including  ini¬ 
tiation,  termination,  propagation,  and  branching  re¬ 
actions.  We  can  define  ignition  by  requiring  that  the 
reacting  system  experience  exponential  growth  in 
both  temperature  and  number  of  chain  carriers,  and 
the  exponential  chain  reaction  must  proceed  for  a 
significant  degree  of  fuel  consumption. 

Chain  Reactions 

Early  analyses  [5,6]  established  the  chain  charac¬ 
ter  of  reaction  of  the  H2-02  system,  and  the  same 
principles  apply  to  hydrocarbon  oxidation.  Chain 
carriers  are  radical  species  such  as  H  and  O  atoms, 
OH,  CH3,  and  H02,  and  it  is  usually  straightforward 
to  identify  reactions  as  initiation,  propagation, 
branching,  and  termination. 

Initiation  reactions  generate  radicals  from  stable 
species,  such  as  the  decomposition  of  propane: 

C3H8  -a  CH3  +  C2H5  (3) 

Chain  propagation  reactions  maintain  the  number  of 
radical  species,  as  in 

C2H6  +  OH  -»  C2Hs  +  H20  (4) 

consuming  OH  and  producing  ethyl  radicals.  Chain 
termination  reduces  the  number  of  radicals,  as  in 
recombination  producing  stable  butane: 

C2HS  +  C2H5  -a  C4H10  (5) 

The  key  to  understanding  ignition  kinetics  is  to 
identify  the  chain  branching  steps  under  the  condi¬ 
tions  being  studied.  In  chain  branching  reactions, 
the  number  of  radicals  increases,  as  in 

CH4  +  O  -4  CH3  +  OH  (6) 

consuming  one  O  atom  and  producing  two  radicals. 
The  most  important  high  temperature  chain  branch¬ 
ing  reaction  consumes  one  II  atom  and  produces 
two  radicals,  O  and  OH: 

H  +  02  -4  O  +  OH  (7) 

These  illustrations  are  encapsulated  into  single  re¬ 
actions.  However,  it  is  common  to  find  situations  in 
which  chain  properties  are  best  attributed  to  a  se¬ 
quence  of  reactions.  For  example,  an  important  re¬ 
action  sequence  is 

C3H8  +  II02  -a  C3H7  +  H202  (8) 

H202  -a  OH  +  OH  (9) 

consuming  one  hydroperoxy  radical  and  producing 
three  radicals,  a  propyl  radical  and  two  hydroxyl  rad¬ 
icals.  In  some  situations,  the  sequence  of  reactions 


can  be  much  longer  before  its  chain  properties  be¬ 
come  evident. 

Finally,  not  all  radical  species  have  the  same  influ¬ 
ence  on  the  overall  reaction  rate,  due  to  subsequent 
reactions  with  different  branching  characteristics. 
For  example,  at  high  temperatures,  where  reaction 
7  is  the  major  chain-branching  step,  production  of 
ethyl  radicals  accelerates  the  overall  rate  of  ignition 
because  C2HS  produces  H  atoms,  which  then  pro¬ 
vide 

C2Hs  -4  C2H4  +  H  (10) 

chain  branching  via  reaction  7,  while  production  of 
methyl  radicals  actually  decreases  or  inhibits  the 
overall  rate  of  ignition  at  high  temperatures  because 
methyl  radicals  recombine,  resulting  in  chain  ter¬ 
mination: 

ch3  +  ch3  -4  C2H6  (11) 

The  specific  reaction  sequences  that  provide  chain 
branching  change  as  the  temperature,  pressure,  and 
reactant  composition  change.  In  discussions  below, 
chemical  kinetic  factors  for  ignition  under  a  variety 
of  conditions  will  be  examined,  and,  for  each  envi¬ 
ronment,  the  essential  task  will  be  to  identify  the 
chemical  kinetic  reaction  sequences  that  lead  to 
chain  branching. 

High  Temperature  Ignition,  Shock  Tubes, 
Detonations,  and  Pulse  Combustors 

At  temperatures  above  about  1200  K,  the  domi¬ 
nant  chain-branching  step  in  hydrocarbon  ignition  is 
reaction  7  [7-12],  with  H  atoms  that  are  produced 
by  thermal  decomposition  of  radicals  such  as  ethyl, 
vinyl,  formyl,  isopropyl,  and  others.  Activation  en¬ 
ergies  of  these  reactions  are  quite  large  (e.g.,  30 
kcal/mol),  and  the  activation  energy  of  reaction  7  is 
also  quite  large  (16.8  kcal/mol),  so  these  reactions 
need  high  temperatures  to  proceed. 

Shock  Tubes 

High  temperature  (T  a  1200  K)  shock  tube  igni¬ 
tion  is  controlled  by  reaction  7.  At  shock  tube  tem¬ 
peratures,  there  is  an  initiation  period  for  radical 
generation,  followed  by  a  period  of  explosive,  ex¬ 
ponential  growth  in  radicals,  dominated  by  reaction 
7,  consuming  fuel  and  increasing  the  temperature. 
An  interesting  series  of  experiments  was  reported  by 
Burcat  et  al.  [13]  for  ignition  of  stoichiometric 
mixtures  of  n-alkanes  in  oxygen/argon.  Fuels  in¬ 
cluded  methane,  ethane,  propane,  /i -butane,  and  n- 
pentane.  Measured  ignition  delay  times  are  sum¬ 
marized  as  symbols  in  Fig.  1,  together  with 
computed  values  using  kinetic  reaction  mechanisms, 
with  the  exception  of  ethane;  the  heavy  dashed  line 
is  a  fit  to  the  experimental  results  for  ethane  ignition, 
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Fig.  1.  Ignition  delay  times  for  C^Cs  alkanes.  Experi¬ 
mental  data  are  from  Burcat  et  al.  [13];  computed  lines  are 
Westbrook  and  Pitz  [10],  The  dark  dashed  line  is  a  fit  to 
the  experimental  results  for  ethane  ignition  from  Ref.  [13], 
and  no  computed  results  for  ethane  ignition  were  carried 
out.  Experimental  points:  open  triangles,  methane;  filled 
circles,  propane;  filled  squares,  n-pentane;  open  squares, 
n-butane. 

and  no  computed  results  for  ethane  are  presented 
in  Fig.  1.  Calculated  results  for  isobutane  are  indi¬ 
cated  as  a  dotted  line  and  are  virtually  indistinguish¬ 
able  from  n-butane  and  n-pentane.  Of  these  five  n- 
alkane  fuels,  methane  ignites  slowest,  ethane  most 
rapidly,  and  the  three  higher  n-alkanes  have  inter¬ 
mediate  reactivity  and  are  very  similar  to  each  other. 
Kinetic  modeling  shows  that  methane  ignites  slowest 
because  methyl  radicals,  resulting  from  H-atom  ab¬ 
straction  from  methane,  lead  to  chain  termination 
through  reaction  11.  Ethane  is  most  reactive  because 
every  ethyl  radical,  resulting  from  H-atom  abstrac¬ 
tion  from  ethane,  produces  an  H  atom  via  reaction 
10.  For  higher  n-alkanes,  two  or  more  alkyl  radicals 
are  produced,  some  of  which  produce  II  atoms  and 
chain  branching  through  reaction  7  and  some  of 
which  produce  methyl  radicals  and  chain  termina¬ 
tion  through  reaction  11.  For  example,  butyl  radicals 
from  n-butane  decompose: 

1- C4H9  -a  C2H4  +  C2H5  (12) 

2- C4H9  -a  C3H6  +  CH3  (13) 

Since  reaction  10  produces  H  atoms  from  C2H5, 
these  two  steps  show  that  n-butane  produces  a  mix¬ 
ture  of  H  atoms  and  CH3  radicals,  and  the  same  is 
true  with  propane  and  n-pentane.  Thus,  methane 
and  ethane  represent  extremes  of  reactivity,  while 
larger  alkane  hydrocarbons  have  intermediate  igni¬ 
tion  properties  due  to  their  mixed  production  of  H 
and  CH3. 


Sensitivity  of  the  high  temperature  mechanism  to 
reaction  7  is  responsible  for  many  practical  phenom¬ 
ena.  Flame  and  detonation  inhibitors  catalyze  re¬ 
moval  of  H  atoms  [14-19],  as  illustrated  here  in  the 
case  of  HBr: 


H  +  HBr  - 

4  H2  +  Br 

(14) 

H  +  Br2  - 

4  HBr  +  Br 

(15) 

Br  +  Br  +  M  - 

4  Br2  +  M 

(16) 

for  a  net  reaction  of 

H  +  H  — >  H2 

This  net  removal  of  I  I  atoms  reduces  the  number  of 
H  atoms  available  for  chain  branching  via  reaction  7 
and  slows  the  overall  rate  of  ignition. 

This  discussion  addresses  the  most  common  shock 
tube  experiments.  However,  some  shock  tube  ex¬ 
periments  are  carried  out  at  lower  temperatures  [20] 
or  at  significantly  elevated  pressures  [21],  where  the 
chain  branching  reaction  sequences  are  different 
from  the  conventional  high  temperature  conditions 
outlined  here;  these  reaction  environments  will  be 
discussed  below. 

Detonations 

High  temperature  ignition  kinetics  and  reaction  7 
play  essential  roles  in  the  propagation  of  gaseous 
detonations,  where  postcompression  detonation 
temperatures  and  pressures  are  comparable  with 
those  in  shock  tubes.  Although  it  is  often  convenient 
to  consider  a  detonation  as  a  planar  phenomenon, 
detonations  have  a  complex,  three-dimensional 
structure  caused  by  interactions  between  transverse 
waves  in  the  reactive  medium  [22,23],  The  reactive 
shock  wave  decays  during  propagation  and  would  fail 
if  new  ignition  spots  were  not  created  by  intersecting 
shock  waves,  or  Mach  stems.  For  this  reason,  deto¬ 
nation  propagation  depends  veiy  strongly  on  ignition 
kinetics.  The  result  of  these  interactions  is  a  cellular 
structure  that  can  be  measured  using  smoked  foils 
in  laboratory  experiments.  The  widths  of  these  det¬ 
onation  cells  are  fundamental  features  relating  to  the 
initiation,  propagation,  and  stability  of  gaseous  det¬ 
onations. 

Because  detonation  structure  is  inherently  three- 
dimensional,  and  the  reaction  zone  is  spatially  thin 
relative  to  cell  sizes,  computer  models  for  detona¬ 
tions  [24]  have  been  unable  to  include  full  chemical 
submodels,  although  recent  models  [25]  have  used 
tabular  look-up  techniques  to  include  finite-rate  ki¬ 
netics.  Advances  in  massively  parallel  computer  ca¬ 
pabilities  are  likely  to  permit  inclusion  of  detailed 
kinetics  in  future  three-dimensional  CFD  models  of 
detonations. 

Given  the  computer  costs  of  a  three-dimensional 
detonation  model,  the  pseudo-one-dimensional 


1566 


INVITED  TOPICAL  REVIEW 


Zeldovich-von  Neumann-Doring  (ZND)  model  [26- 
28]  has  been  used  to  relate  computed  ignition  delay 
times  to  characteristic  spatial  and  energy  scales  in  a 
detonation.  This  approach  calculates  ignition  delay 
times  behind  a  Chapman-Jouguet  (CJ)  shock  wave 
to  define  a  characteristic  length  scale  in  the  reacting 
gas  mixture  [29-37],  which  is  then  proportional  to 
the  cell  width,  critical  tube  diameter,  minimum  ini¬ 
tiation  energy,  and  other  detonation  characteristics. 

An  exciting  current  development  is  extension  of 
kinetic  modeling  to  study  detonations  of  condensed 
phase  explosives  [38-41],  such  as  ammonium  per¬ 
chlorate,  RDX,  or  HMX.  These  models  require 
treatment  of  three-dimensional  fluid  mechanics, 
multiphase  phenomena,  and  detailed  chemical  ki¬ 
netics,  and  massively  parallel  computing  is  an  essen¬ 
tial  computational  tool.  Ignition  kinetics  in  such  an 
environment  are  complex,  and  a  great  deal  of  new 
research  is  needed. 


Pulse  Combustion 

Pulse  combustion  has  been  known  for  centuries, 
but  only  recently  have  its  physical  and  chemical  prin¬ 
ciples  become  understood.  Ignition  and  heating  in 
the  pulse  combustor  force  exhaust  gases  out  through 
an  exhaust  pipe;  the  resulting  lower  pressure  in  the 
combustion  chamber  then  draws  in  fresh  fuel  and 
air.  This  fresh  mixture  combines  with  hot  residual 
products  of  the  previous  burn,  and  the  resulting 
fresh  reactants  and  residual  gases  ignite,  starting  a 
new  cycle.  Because  the  average  temperature  is  rela¬ 
tively  low,  production  of  NOr  is  also  low. 

Keller  et  al.  [42-45]  examined  mixing,  kinetics  of 
ignition,  and  resonant  acoustic  wave  propagation  in 
the  pulse  combustor,  showing  that  ignition  must  be 
timed  to  occur  in  phase  with  resonant  pressure 
waves  in  the  combustion  chamber. 

The  kinetic  submodel  for  this  system  mixes  fresh 
fuel  and  air,  initially  at  room  temperature,  with  hot 
residual  gases  from  the  preceding  cycle,  steadily  in¬ 
creasing  the  temperature  of  the  fresh  mixture  as  well 
as  diluting  it  with  the  residual  products.  The  overall 
reaction  rate  remains  negligible  until  the  mixture 
reaches  temperatures  of  about  1200  K,  where  igni¬ 
tion  is  controlled  by  reaction  7.  The  kinetic  model 
for  this  system  predicts  the  influence  of  additives  and 
variability  in  natural  gas  composition  [46]  to  manip¬ 
ulate  the  kinetics  of  ignition  in  a  pulse  combustor. 
This  is  a  veiy  rich  system  in  terms  of  the  parameters 
that  can  be  adjusted  to  optimize  the  phasing  be¬ 
tween  kinetic  ignition  and  the  resonant  pressure 
wave  in  the  combustion  chamber.  In  addition  to  the 
kinetic  parameters,  the  combustion  chamber  acous¬ 
tic  timescales  can  be  modified,  and  the  mixing  rate 
can  be  varied  by  adjusting  intake  ports  for  fuel  and 
air. 


Intermediate  Temperature  Ignition,  Rapid 
Compression  Machines,  Engine  Knock,  Diesel 
Ignition,  and  Homogeneous  Charge, 
Compression  Ignition  (HCCI) 

At  temperatures  above  about  850  K  but  below 
1200  K,  reaction  7  is  too  slow  to  provide  sufficient 
branching  rates  for  ignition,  and  a  different  reaction 
path  dominates.  Key  reactions  include  the  following: 


+  O2  +  M  — 

A  HOa  +  M 

(17) 

RH  +  H02  - 

->  R  +  H2O2 

(18) 

H202  +  M  — 

->  OH  +  OH  +  M 

(19) 

where  RH  is  an  alkane,  R  is  an  alkyl  radical,  and  M 
is  a  third  body.  Collectively,  reactions  17-19  con¬ 
sume  one  H-atom  radical  and  produce  two  OH  rad¬ 
icals,  providing  chain  branching. 

Kinetic  studies  in  flow  reactors  and  stirred  reactors 
[47-49]  at  temperatures  of  about  1000  K  show  that 
reactions  17-19  control  combustion  rates.  In  this 
range,  reaction  19  is  rapid,  so  H202  decomposes  as 
rapidly  as  it  is  formed  and  H202  concentrations  do 
not  build  up  to  appreciable  levels.  However,  in  many 
practical  systems,  the  dominant  feature  of  ignition  is 
the  accumulation  of  H202  until  a  temperature  is 
reached  where  it  decomposes  to  cause  ignition.  This 
is  true  in  the  laboratory  rapid  compression  machine, 
and  it  is  also  the  central  kinetic  feature  in  engine 
knock  in  spark  ignition  engines,  in  ignition  in  liquid- 
fueled  diesel  engines,  and  in  the  operation  of  ho¬ 
mogeneous  charge,  compression  ignition  (HCCI) 
engines.  In  each  of  these  systems,  H202  forms  at 
lower  temperatures  and  remains  relatively  inert,  un¬ 
til  increasing  temperature  from  compression  and  ex¬ 
othermic  reactions  reaches  a  level  where  it  decom¬ 
poses  rapidly  via  reaction  19.  Each  of  these  examples 
will  be  discussed  below. 

Rapid  Compression  Machine 

Rapid  compression  machines  (RCMs)  have  been 
used  for  many  years  [50-52]  and  provide  a  fertile 
environment  for  the  study  of  hydrocarbon  oxidation. 
The  apparatus  provides  a  combustion  chamber  with 
a  single  piston  stroke,  leaving  the  compressed  gases 
at  temperatures  from  550  to  950  K,  depending  on 
the  compression  ratio  of  the  RCM  and  specific  heats 
of  reactant  gases.  Compressed  gases  can  undergo  a 
single-stage  or  two-stage  ignition,  or  the  reactants 
may  not  ignite  at  all  while  heat  transfer  cools  the 
mixture  to  room  temperature. 

Typical  results  in  Fig.  2  are  taken  from  numerical 
simulation  of  ignition  of  a  stoichiometric  mixture  of 
neopentane  and  oxygen  in  N2  and  Ar  in  proportions 
of  CsH12/02/N2/Ar  =  0.0256/0.2048/0.38/0.39  at 
a  postcompression  temperature  of 757  K  [53],  In  this 
example,  the  temperature  history  shows  a  well-de¬ 
fined  two-stage  ignition,  the  first  stage  occurring  af¬ 
ter  a  delay  of  5  ms,  reaching  a  new  temperature  of 
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Fig.  2.  Computed  temperature  and  concentrations  of 
neopentane,  H202,  and  OH  from  Ribaucour  et  al.  [53], 


about  880  K.  The  temperature  then  remains  rela¬ 
tively  constant,  rising  slowly  for  the  next  6-8  ms, 
then  rising  more  rapidly  until  a  total  elapsed  time  of 
14  ms,  when  the  second-stage  ignition  is  observed 
and  the  temperature  rises  rapidly,  consuming  all  of 
the  fuel.  In  other  cases,  the  first  stage  can  occur  dur¬ 
ing  the  compression  stroke  [54,55]  and  is  completed 
before  the  end  of  the  compression,  while  in  other 
cases  there  is  no  discernible  first-stage  ignition. 

The  first-stage  ignition  is  a  result  of  rather  complex 
kinetic  features  of  low  temperature  hydrocarbon  ox¬ 
idation  [56],  This  regime  has  considerable  impor¬ 
tance  in  theoretical  and  practical  systems,  and  we 
wall  return  to  describe  it  in  more  detail  below. 

Figure  2  plots  the  concentrations  of  fuel,  H202, 
and  OH.  During  the  time  preceding  the  second- 
stage  ignition,  H202  concentration  increases  stead¬ 
ily.  Until  the  temperature  approaches  1000  K,  de¬ 
composition  of  H202  is  much  slower  than  its 
production,  and  the  small  amounts  of  OH  produced 
are  consumed  by  reactions  with  the  fuel  [57],  When 
a  temperature  close  to  1000  K  is  reached,  reaction 
19  accelerates,  H202  disappears,  and  OH  concen¬ 
tration  increases  rapidly.  The  increased  OH  concen¬ 
trations  rapidly  consume  any  remaining  fuel,  fol¬ 
lowed  by  a  rapid  increase  in  temperature.  Thus, 
decomposition  of  H202  and  consumption  of  fuel  re¬ 
sult  in  ignition. 

There  is  a  critical  temperature  for  ignition  through 
H202  decomposition  that  is  a  function  of  the  pres¬ 
sure  of  the  reactive  system.  This  is  seen  by  examining 
the  rate  of  H202  decomposition 

dlHft  h]  =  -  [H202]  [M ]k  (20) 

where  M  is  the  total  molar  concentration.  This  equa¬ 
tion  can  be  rearranged  into  the  form  [H202]/ 
(d[H202]/c?f);  a  characteristic  decomposition  time  is 

r  =  [H202]/(d[H202]/dt)  =  1/(*[M])  (21) 

With  the  rate  constant  for  reaction  19,  =  1.2  X 

1017  *  exp(  —45,500/117’),  r  becomes 


r  =  8.3  X  10“18  * 

exp(  +  22,750/T)  *  [M]_1  (22) 

As  the  temperature  increases,  r  becomes  smaller, 
and  when  it  is  “short”  compared  with  the  residence 
time,  ignition  occurs.  Characteristic  decomposition 
time  also  decreases  with  increasing  total  concentra¬ 
tion  [M]  or  with  increasing  pressure  at  constant  tem¬ 
perature.  Therefore,  as  pressure  increases,  the  criti¬ 
cal  temperature  for  ignition  decreases  gradually. 

Total  postcompression  concentrations  in  a  repre¬ 
sentative  RCM  experiment  are  1  X  10-4  mol/cm3, 
so  r  is  approximately  7.8  ms  at  T  —  900  K,  640  /is 
at  T  =  1000  K,  and  80 /is  at  T  =  1100  K.  With  total 
reaction  times  of  tens  of  milliseconds,  rapid  decom¬ 
position  is  observed  at  about  950-1000  K,  with  char¬ 
acteristic  times  for  ignition  less  than  a  millisecond. 
In  studies  of  ignition  at  the  higher  pressures  of  in¬ 
ternal  combustion  engines,  including  spark  ignition, 
diesel,  and  HCCI  engines,  the  same  calculation  of 
equations  20-22  shows  that  this  simple  mechanism 
predicts  ignition  at  temperatures  between  900  and 
950  K  [58],  Analysis  of  reaction  19  by  Griffiths  and 
Barnard  [1]  at  its  high  pressure  limit  shows  similar 
rates  of  decomposition  at  somewhat  lower  tempera¬ 
tures,  consistent  with  the  above  trends  as  pressure 
increases. 

The  key  is  that  H202  decomposes  rapidly  at  tem¬ 
peratures  of  900-1000  K,  so  this  temperature  can  be 
considered  as  a  critical  value  in  systems  in  which  the 
temperature  approaches  this  value  from  lower  tem¬ 
peratures.  H202  is  produced  at  lower  temperatures 
by  low  and  intermediate  temperature  kinetic  path¬ 
ways,  until  the  decomposition  temperature  is 
reached,  when  the  system  suddenly  produces  large 
numbers  of  OH  radicals  and  rapidly  ignites.  The 
most  important  variable  is  the  time  at  which  they 
reach  this  critical  temperature.  Anything  that  wall  ac¬ 
celerate  reaching  this  critical  temperature  advances 
its  ignition,  while  delaying  it  inhibits  ignition.  The 
importance  of  the  first-stage,  low  temperature  oxi¬ 
dation  period  is  that  it  provides  heat  release  early  in 
the  reaction  history,  so  the  reactive  mixture  arrives 
at  the  decomposition  temperature  earlier  than 
would  have  occurred  without  that  early  heat  release. 

Low  Temperature  Chemical  Kinetics 

Low  temperature  hydrocarbon  oxidation  has  been 
studied  extensively  [56] .  F ollowing  abstraction  of  H 
atoms  from  the  fuel,  at  high  temperatures  the  alkyl 
radical  R  decomposes,  producing  olefin  and  smaller 
alkyl  radicals,  and  reaction  7  is  the  dominant  chair¬ 
branching  reaction,  as  discussed  above.  At  lower 
temperatures,  02  adds  to  the  alkyl: 

R  +  02  +  M  -A  ROz  +  M  (23) 
The  equilibrium  constant  for  reaction  23  is  strongly 
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Fig.  3.  Computed  temperatures  in  ignition  of  the  three 
isomers  of  pentane,  from  Ribaucour  et  al.  [53]. 

temperature  dependent  and  is  strongly  in  favor  of 
ROa  at  low  temperature,  shifting  toward  R  +  Oz  as 
the  temperature  increases.  R02  radicals  isomerize  to 
produce  a  QOOH  radical  species 

ROz  QOOH  (24) 

Isomerization  depends  sensitively  on  the  size  and 
structure  of  the  original  fuel  molecule  and  the  site 
in  that  fuel  where  the  02  group  is  located  [56,58]. 
The  R02  will  be  produced  with  excess  energy,  so  it 
is  important  to  understand  the  rate  of  collisional  sta¬ 
bilization  and  whether  subsequent  reactions  occur 
from  the  activated  or  ground  state  of  R02.  The  sim¬ 
plest  example  of  such  a  system  is  the  ethyl  radical, 
with  R02  being  the  ethylperoxy  radical  C2H502  and 
QOOH  being  the  hydroperoxyethyl  radical 
C2H4OOH.  This  system  has  been  studied  exten¬ 
sively  [59-63],  and  even  this  small  system  is  not  fully 
understood.  Bozzelli  and  Pitz  [64]  and  Koert  et  al. 
[65]  examined  the  propyl  +  02  system,  but  little 
experimental  or  theoretical  work  has  been  done  for 
larger  systems. 

Several  factors  influence  the  rate  of  internal  H- 
atom  transfer  or  radical  isomerization.  The  species 
forms  a  ringlike  transition  state  where  the  terminal 
O  atom  approaches  an  extractable  H  atom  within  the 
R02  species.  The  number  of  atoms  in  this  transition 
state  ring  influences  the  rate  of  reaction,  with  six- 
and  seven-membered  rings  having  the  most  rapid 
rate.  The  type  of  C-H  bond  broken  has  a  strong  in¬ 
fluence,  with  primary  C-H  bonds  being  strongest, 
secondary  C-H  bonds  next,  and  tertiary  C-H  bonds 
being  weakest.  The  rate  and  number  of  possible  R02 
isomerization  reactions  increase  with  the  fuel  mol¬ 
ecule  size  and  are  fastest  in  long,  linear  alkane  fuel 
molecules,  with  a  large  number  of  six-  and  seven- 
membered  transition  state  rings  and  a  high  percent¬ 
age  of  easily  abstracted  secondary  C-H  bonds,  and 
are  slowest  in  highly  branched  fuel  molecules,  with 
fewer  low-energy  transition  state  rings  and  large  per¬ 
centages  of  difficult-to-abstract  primary  C-H  bonds 


QOOH  species  react  via  several  alternative  paths 
which  are  formally  chain  propagation  steps, 

QOOH  ->  Q  +  H02  (25) 

QOOH  ->  QO  +  OH  (26) 

beginning  with  one  alkyl  radical,  and  produce  one 
H02  or  OH  radical,  in  addition  to  a  stable  olefin  or 
cyclic  ether.  However,  it  is  also  possible  for  a  second 
02  molecule  to  add  to  QOOH,  creating  a  radical 
02QOOH,  which  can  then  isomerize  further 

QOOH  +  02  02QOOH  (27) 

again  with  rates  dependent  on  the  structure  and  size 
of  the  original  fuel  molecule,  as  observed  for  ROa 
isomerization.  The  isomerized  product  decomposes 
into  a  relatively  stable  ketohydroperoxide  species 
and  one  OH  radical.  The  ketohydroperoxide  species 
then  has  its  own  temperature  for  decomposition  at 
about  800  K,  somewhat  lower  than  that  of  H202. 
Upon  reaching  this  temperature,  the  ketohydrope¬ 
roxide  decomposes  into  several  pieces,  at  least  two 
of  which  are  radicals.  Thus,  it  is  not  until  this  final 
ketohydroperoxide  decomposition  step  that  chain 
branching  is  finally  achieved  in  the  low  temperature 
oxidation  regime.  Since  at  least  three  of  the  ultimate 
products  of  this  reaction  sequence  are  radicals,  chain 
branching  is  quite  strong  once  the  ketohydroperox¬ 
ide  decomposes.  Production  of  metastable  inter¬ 
mediates,  ketohydroperoxide  in  this  case  and  H202 
earlier,  followed  later  by  the  higher  temperature  de¬ 
composition  of  the  metastable  intermediate  to  finally 
provide  chain  branching  are  examples  of  degenerate 
chain  branching.  Detailed  low  temperature  reaction 
mechanisms  have  been  developed  [53-55,58,66,67] 
for  a  variety  of  hydrocarbon  fuels. 

This  highly  branched,  low  temperature  oxidation 
phase  continues  until  the  temperature  has  increased 
enough  that  the  equilibria  in  the  molecular  oxygen 
addition  reactions  23  and  27  begin  to  shift  toward 
dissociation.  Because  activation  energies  for  disso¬ 
ciation  are  large  (i.e.,  ~35  kcal/mol),  the  reactions 
shift  rapidly  toward  dissociation  as  temperature  in¬ 
creases.  This  shift  with  hydrocarbon  fuels  occurs  at 
temperatures  near  850  K.  The  temperature  above 
which  these  equilibria  favor  dissociation  has  been 
termed  the  “ceiling  temperature”  by  Benson  [68,69], 

The  effect  of  fuel  molecular  structure  on  both  the 
first-  and  second-stage  ignition  is  illustrated  in  Fig. 
3,  from  Ribaucour  et  al.  [53],  The  three  isomers  of 
pentane  are  compared  at  the  same  postcompression 
temperature  of  757  K.  All  three  mixtures  ignite  at 
about  950  K.  n-pentane  ignites  first  since  its  first 
stage  occurs  first  and  provides  the  largest  tempera¬ 
ture  increase,  neopentane  (2,2-dimethyl  propane) 
ignites  next,  and  isopentane  (2-methyl  butane)  is 
last,  all  because  the  time  of  occurrence  and  tem¬ 
perature  increase  of  the  first-stage  ignition  vaiy  in 
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Time  after  TDC  [ms] 

Fig.  4.  Computed  temperature  and  selected  species 
concentrations  from  rapid  compression  machine  ignition 
of  neopentane.  Same  initial  conditions  as  Fig.  2,  but  with 
10  ppm  ozone  included. 

that  order.  All  ignite  when  they  reach  the  same  criti¬ 
cal  temperature  for  II202  decomposition,  and  the 
differences  between  the  isomers  are  the  times  when 
the  fuel  mixtures  reach  the  critical  temperature. 

Total  heat  release  in  this  low  temperature  regime 
is  quite  modest.  In  Fig.  3,  this  varies  from  50  K  for 
2-methyl  butane  to  150  K  for  n -pentane.  In  engine 
studies,  at  top  dead  center  (TDC)  this  difference 
makes  a  significant  impact  on  ignition  properties  of 
the  combustion  system. 

An  interesting  variation  of  this  analysis  simulates 
an  RCM  experiment  in  which  an  additive  has  been 
included  in  the  neopentane  fuel.  This  additive, 
ozone,  has  an  initial  concentration  of  10  ppm  and 
decomposes  early  in  the  experimental  cycle,  since 
the  activation  energy  for  ozone  decomposition  is  24 
kcal/mol,  much  less  than  the  45.5  kcal/mol  for  H202 
decomposition.  Repeating  the  analysis  of  equations 
20-22,  the  critical  decomposition  temperature  for 
ozone  at  present  conditions  is  about  600  K.  Com¬ 
puted  results  of  this  example  are  shown  in  Fig.  4. 

Decomposition  of  ozone  at  5  ms  before  TDC  (at 
~600  K)  produces  O  atoms.  Each  O  atom  consumes 
two  fuel  molecules 

03  +  M^02  +  0  +  M  (28) 

CsH12  +  O  ->  C5Hu  +  OH  (29) 

C5H12  +  OH  C5Hu  +  II20  (30) 

Pentyl  radicals  then  proceed  via  reactions  23-27  and 
produce  the  metastable  C5  ketohydroperoxide  spe¬ 
cies,  the  production  of  which  is  shown  in  Fig.  4,  im¬ 
mediately  following  the  decomposition  of  ozone. 

Ketohydroperoxide  decomposition  has  an  activa¬ 
tion  energy  of  about  42  kcal/ mol,  so  it  decomposes 
at  about  800  K,  as  seen  in  Fig.  4.  Production  of  OH 
from  decomposition  produces  further  water  and 
heat  release.  Decomposition  of  H202  and  the  keto¬ 
hydroperoxide  and  the  resulting  heat  release  bring 


the  reactive  mixture  to  the  first  ignition  stage  before 
the  conclusion  of  the  compression  stroke.  Fig.  4 
shows  a  slight  inflection  point  at  TDC,  indicating 
that  the  low-temperature  ignition  is  well  under  way 
as  the  compression  stroke  ends.  This  phase  stops  at 
about  850  K,  the  same  temperature  at  which  the  first 
stage  ended  in  the  original  mixture,  without  ozone 
(Figs.  2  and  3).  The  kinetics  of  the  first  stage  have 
not  been  altered  by  the  additive,  but  the  time  at 
which  it  begins  was  advanced  by  adding  ozone, 
which  stimulated  early  heat  release. 

The  major  intermediate  being  produced  is  still 
H202,  and  it  still  decomposes  at  950-1000  K,  pro¬ 
ducing  the  real  ignition.  The  ozone  additive  ad¬ 
vances  the  time  of  ignition  from  about  14  ms  after 
TDC  to  about  9  ms  after  TDC.  The  second-stage 
ignition  is  the  same  in  both  cases,  but  the  addition 
of  ozone  makes  that  mixture  reach  the  ignition  tem¬ 
perature  at  an  earlier  time.  Finally,  note  the  peak  in 
O-atom  concentration  that  occurs  at  each  of  the 
three  stages  of  ignition,  at  times  of  —  6,  +1,  and  +  9 
ms  after  TDC. 

Engine  Knock 

In  spark  ignition  engines,  thermodynamic  com¬ 
bustion  efficiency  increases  with  engine  compres¬ 
sion  ratio.  However,  an  increased  compression  ratio 
eventually  results  in  engine  knock,  limiting  engine 
efficiency  in  practical  engines.  Numerical  studies 
have  used  reduced  [70-74]  and  detailed  [66,74-78] 
kinetic  models  to  understand  chemical  factors  re¬ 
sponsible  for  knock. 

In  a  spark  ignition  engine,  a  flame  propagates 
through  a  combustion  chamber,  starting  at  the  point 
of  spark  ignition.  Some  reactant  gases  in  the  com¬ 
bustion  chamber  will  naturally  be  the  last  to  be  con¬ 
sumed  by  the  flame,  termed  the  end  gases.  End  gas 
conditions  are  determined  by  piston  motion  and 
combustion  in  the  engine  chamber.  As  piston  motion 
and  flame  propagation  proceed,  end  gases  see  in¬ 
creasing  levels  of  pressure  and  temperature  and  re¬ 
act  accordingly,  if  end  gases  reach  the  point  of  ig¬ 
nition  before  being  consumed  by  the  flame, 
knocking  behavior  will  be  observed.  Mixtures  that 
react  more  rapidly  are  more  susceptible  to  knock, 
while  mixtures  that  ignite  more  slowly  resist  knock. 
Increasing  the  engine  compression  ratio  increases 
the  rate  of  autoignition  while  having  little  effect  on 
flame  propagation,  so  this  increases  the  potential  for 
knocking  behavior.  Thus,  engine  knock  is  essentially 
an  ignition  problem. 

Computational  modeling  of  autoignition  at  engine 
conditions  shows  the  same  features  as  described  for 
the  rapid  compression  machine.  This  can  be  illus¬ 
trated  by  examining  the  way  that  fuel  composition 
and  molecular  size  and  structure  influence  autoig¬ 
nition  behavior  and  tendency  to  knock  by  simulating 
the  critical  compression  ratio.  Experiments  were 
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600  rpm,  T=403  K,  MAP=81.6  kPa 


Fig.  5.  Correlation  of  computed  critical  compression  ra¬ 
tio  with  research  octane  numbers.  Solid  circles  are  primary 
reference  fuel  mixtures,  other  symbols  are  isomers  of  hex¬ 
ane  and  pentane,  and  stars  are  experimentally  measured 
values  [66], 
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Fig.  6.  Schematic  model  of  diesel  spray  combustion, 
based  on  the  work  of  Dec  [83] .  The  ignition  region  is  iden¬ 
tified,  and  the  temperatures  of  the  different  regions  in  the 
combusting  spray  are  noted. 

carried  out  in  a  Cooperative  Fuels  Research  engine 
at  a  low  engine  speed  of  600  rpm  and  intake  mani¬ 
fold  temperature  of  403  K.  These  are  essentially  the 
same  conditions  as  those  in  standardized  tests  for 
research  octane  number  (RON).  For  stoichiometric 
fuel/air  mixtures,  the  engine  compression  ratio  is 
steadily  increased  until  autoignition  is  observed.  A 
series  of  such  experiments  were  carried  out  by  Lep¬ 
pard  [79]  for  a  range  of  alkane  and  primary  reference 
fuel  mixtures,  and  a  kinetic  model  was  used  to  sim¬ 
ulate  those  experiments  [66,72],  Computed  and  ex¬ 
perimental  critical  compression  ratio  results  are 
plotted  in  Fig.  5,  together  with  the  measured  RONs 
of  each  fuel.  A  rather  smooth  curve  results,  although 
it  is  clearly  not  a  straight  line,  indicating  that  octane 
rating  is  a  rather  nonlinear  scale. 


A  numerical  experiment  can  be  carried  out  to  use 
this  information  to  construct  a  fuel  mixture  with  an 
arbitrary  octane  rating.  For  example,  to  create  a  fuel 
with  a  RON  of  90.8,  a  motor  octane  number  (MON) 
of  83.9,  and  an  overall  octane  rating  (RON  + 
MON)/2  of  89,  to  compare  with  ordinary  gasoline, 
a  mixture  of  two  isomers  of  hexane  and  two  isomers 
of  pentane  was  defined  in  terms  of  their  individual 
RON  and  MON  values.  These  four  fuel  components 
are  indicated  by  the  arrows  in  Fig.  5.  When  tested 
numerically,  this  mixture  had  a  critical  compression 
ratio  corresponding  to  a  RON  of  90.4,  veiy  close  to 
the  blending  value,  as  shown  in  Fig.  5. 

A  critical  feature  of  these  computed  histories  is 
that  they  reach  peak  compression  temperatures  that 
are  all  very  nearly  the  same.  Considerable  differ¬ 
ences  in  octane  number  result  in  only  a  few  degrees 
of  temperature  at  TDC,  resulting  in  a  considerable 
difference  in  time  of  ignition  and  demonstrating  that 
the  octane  number  scale  is  highly  nonlinear. 

Just  as  observed  in  the  RCM  results,  these  com¬ 
puted  results  indicate  that  ignition  in  each  of  these 
cases  occurs  when  H202  decomposes.  Differences 
in  octane  number  are  reflected  in  small  differences 
in  cool  flame  heat  release,  with  greater  amounts  of 
low  temperature  heat  release  and  higher  quantities 
of  H202  correlating  with  earlier  ignition  and  lower 
octane  values.  Similar  conclusions  were  obtained  by 
Rlin-Simiand  et  al.  concerning  the  central  role  of 
FLOq  decomposition  and  its  role  in  autoignition 
[80], 

The  effects  of  various  antiknock  compounds  can 
be  understood  in  this  same  framework.  Additives  in¬ 
cluding  tetraethyl  lead,  now  no  longer  used,  provide 
kinetic  sinks  for  H02  [58,81,82],  greatly  reduce  the 
production  of  H202,  and  interfere  with  the  inter¬ 
mediate  temperature  ignition  process.  Other  addi¬ 
tives,  including  methyl  tert -butyl  ether,  act  by  inter¬ 
fering  with  low  temperature  oxidation  [18,58], 
reducing  the  amount  of  low  temperature  heat  re¬ 
lease  and  retarding  the  time  at  which  the  end  gas 
reaches  the  H202  decomposition  temperature. 

Diesel  Ignition 

Diesel  engines  have  existed  for  many  years,  but 
until  recently  many  of  the  basic  physical  and  chem¬ 
ical  principles  of  diesel  combustion  had  not  been 
well  understood.  Recently,  in  a  series  of  insightful 
laser  diagnostic  studies,  Dec  provided  a  coherent, 
self-consistent  picture  of  diesel  combustion  [83],  His 
results,  summarized  in  Fig.  6,  show  that  the  fuel  jet 
vaporizes  rapidly  and  mixes  with  hot,  compressed  air. 
The  air  steadily  reduces  the  fuel/ air  equivalence  ra¬ 
tio  at  the  same  time  that  it  is  increasing  the  mixture 
temperature,  and  the  mixture  begins  to  react.  In 
Dec’s  observations,  this  mixture  eventually  ignites 
while  the  equivalence  ratio  is  still  quite  high  (<j>  ~  4). 
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Fig.  7.  Experimentally  observed  variation  in  Boseh 
smoke  number  with  fuel  oxygenate  content,  from  Miya¬ 
moto  et  al.  [90],  and  variation  in  computed  soot  precursor 
concentrations  with  fuel  oxygenate  content,  from  Flynn  et 
al.  [84], 


Soot  production  was  observed  to  proceed  immedi¬ 
ately  from  the  products  of  this  rich,  premixed  igni¬ 
tion,  and  soot  was  consumed  in  a  diffusion  flame  at 
the  periphery  of  the  downstream  cloud. 

This  large-scale  ignition  problem  has  been  ana¬ 
lyzed  using  detailed  chemical  kinetics  [84],  The  pre¬ 
mixed  region  begins  to  react  when  the  local  equiv¬ 
alence  ratio  reaches  (f>  10,  although  the  rate  of 

reaction  is  initially  quite  slow.  As  air  mixes  and  the 
mixture  temperature  increases,  the  rate  of  reaction 
increases.  Rapid  reaction  begins  when  the  mixture 
temperature  reaches  about  700  K.  As  noted  above 
for  other  classes  of  applications,  this  low  tempera¬ 
ture  reaction  produces  H202  and  produces  a  modest 
amount  of  heat  release  and  temperature  increase. 
Eventually,  the  mixture  reaches  the  temperature  at 
which  H202  decomposes  and  ignition  is  observed. 
Thus,  the  kinetic  mechanism  of  diesel  ignition  is 
identical  to  that  of  the  RCM  and  that  of  engine 
knock  in  spark  ignition  engines.  The  major  differ¬ 
ences  occur  because  diesel  ignition  takes  place  un¬ 
der  very  fuel-rich  conditions. 

Diesel  ignition  improvers  [85]  are  species  such  as 
ethyl-hexyl  nitrate  that  decompose  at  temperatures 
much  lower  than  the  ignition  temperature  provided 
by  H202  decomposition.  Radicals  produced  by  de¬ 
composition  of  the  additive  consume  some  fuel  and 


release  some  heat,  raising  the  temperature  of  the 
premixed  gases  and  getting  them  closer  to  the  H202 
decomposition  temperature.  Other  additives  that 
decompose  at  lower  temperatures  and  provide  rad¬ 
icals  would  also  be  effective  diesel  ignition  improv¬ 
ers  or  cetane  improvers.  The  example  of  ozone  dis¬ 
cussed  earlier  for  the  RCM  would  enhance  diesel 
ignition  in  this  same  manner. 

An  especially  interesting  additional  feature  of  this 
premixed  ignition  is  the  observation  by  Dec  that  the 
products  of  the  rich  premixed  ignition  immediately 
begin  the  process  of  soot  precursor  and  soot  pro¬ 
duction.  In  computed  models  of  this  rich  premixed 
ignition  [84],  the  products  of  the  ignition  are  the 
same  species  that  have  been  shown  [86-89]  to  pref¬ 
erentially  produce  small  aromatic  ring  species  such 
as  benzene,  toluene,  and  naphthalene.  Small  aro¬ 
matic  species  then  react  to  produce  larger  polycyclic 
aromatic  species  and  eventually  soot.  It  was  shown 
experimentally  [90]  that  addition  of  oxygenated  spe¬ 
cies  to  the  fuel  can  reduce  soot  emissions  in  diesel 
combustion,  and  this  trend  was  reproduced  by  ki¬ 
netic  modeling  [84]  by  correlating  the  total  product 
concentrations  of  aromatic  formation  precursors 
with  the  amounts  of  these  simple  species  that  survive 
the  premixed  rich  ignition  stage.  Thus,  soot  produc¬ 
tion  is  directly  related  to  ignition  kinetics  in  the  die¬ 
sel  engine. 

Figure  7  shows  concentrations  of  soot  produced 
in  diesel  engines  with  the  relative  amounts  of  oxy¬ 
genates  in  the  fuel  from  the  experimental  work  of 
Miyamoto  et  al.  [90],  Also  shown  is  the  total  con¬ 
centration  of  ethene,  acetylene,  and  propargyl  radi¬ 
cal,  summed  together,  to  roughly  produce  a  curve 
similar  to  that  derived  experimentally  by  Miyamoto 
et  al.  These  modeling  calculations  [84]  were  carried 
out  for  a  variety  of  oxygenated  species  (i.e.,  metha¬ 
nol,  ethanol,  dimethyl  ether,  dimethoxy  methane, 
and  methyl  butanoate)  added  to  n-heptane;  n-hep- 
tane  is  a  realistic  surrogate  diesel  fuel  with  a  cetane 
number  of  56.  This  work  suggests  strongly  that  rich 
premixed  ignition  does  indeed  produce  the  seeds  of 
soot  formation  and  that  lowering  the  postignition 
concentrations  of  these  elementary  species  effec¬ 
tively  reduces  soot  production.  The  sequence  of 
fuel/air  mixing,  rich  premixed  ignition,  and  produc¬ 
tion  of  soot  precursors,  followed  eventually  by  soot 
burnout,  provides  a  framework  for  diesel  combus¬ 
tion  in  which  each  of  the  detailed  steps  leads  very 
naturally  and  continuously  to  the  next  step. 

Homogeneous  Charge,  Compression  Ignition 

Combustion  in  reciprocating  engines  under  HCCI 
conditions  is  an  intriguing  technology  that  may  offer 
the  opportunity  to  eliminate  N Or  emissions  from  en¬ 
gine  combustion.  Particulate  emissions  are  also  ob¬ 
served  to  be  very  low.  Disadvantages  are  high  un¬ 
burned  hydrocarbon  and  CO  emissions,  along  with 
high  peak  pressures  and  high  rates  of  heat  release. 
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Fig.  8.  Comparison  between  experimental  pressure 
traces  [101]  and  calculated  values,  using  both  a  single-zone 
and  a  10-zone  modeling  treatment  [98],  for  three  natural 
gas  experiments  at  different  boost  pressures. 

In  an  HCCI  engine,  a  premixed,  veiy  lean  mixture 
of  fuel  and  air  is  drawn  into  the  engine  chamber  and 
compressed  by  a  piston.  Near  TDC,  the  majority  of 
the  charge  in  the  engine  ignites  homogeneously.  The 
time  of  ignition  can  be  varied  by  changing  the  com¬ 
pression  ratio,  fuel/air  equivalence  ratio,  time  of  in¬ 
take  valve  opening,  and  intake  manifold  tempera¬ 
ture. 

IICCI  was  identified  as  a  distinct  combustion  phe¬ 
nomenon  about  20  years  ago.  Initial  studies  [91-93] 
recognized  the  basic  characteristics  of  HCCI;  igni¬ 
tion  occurs  at  many  points  simultaneously,  with  no 
flame  propagation.  It  is  controlled  primarily  by 
chemical  ignition  kinetics,  with  little  influence  from 
effects  such  as  turbulence  or  mixing  that  play  such 
a  large  role  in  other  engine  combustion  problems. 
Simplified  kinetics  models  have  been  used  to  analyze 
HCCI  combustion,  coupled  usually  to  multidimen¬ 
sional  CFD  models  [94-98],  with  limited  success. 


Recently,  detailed  kinetic  modeling  was  used  to 
examine  HCCI  combustion  [98-100],  For  a  number 
of  practical  fuels,  including  natural  gas,  propane,  and 
others,  model  calculations  have  reproduced  the  on¬ 
set  of  ignition  in  good  agreement  with  experimental 
results  carried  out  at  the  Lund  Institute  of  Technol¬ 
ogy  [101,102],  using  a  spatially  homogeneous  model. 
The  product  temperature  is  kept  low  by  operating 
at  very  lean  conditions,  so  temperatures  are  too  low 
to  produce  significant  NO,,  emissions.  However,  a 
one-zone,  homogeneous  treatment  produces  much 
too  rapid  a  combustion  duration,  as  shown  in  Fig.  8, 
and  cannot  predict  CO  and  hydrocarbon  emissions. 

HCCI  engine  wall  temperatures  are  low,  so  there 
is  an  extensive  engine  chamber  thermal  boundary 
layer,  and  not  all  the  boundary  layer  fuel  burns  to¬ 
gether  with  the  initial  bulk  ignition.  In  addition, 
there  is  unbumed  fuel  in  the  piston  ring  and  other 
crevice  volumes  that  does  not  burn  immediately. 
Since  the  peak  bulk  temperatures  are  low  and  de¬ 
crease  further  after  TDC,  most  fuel  in  the  boundary 
layer  and  crevices  cannot  diffuse  out  into  the  bulk 
gas  and  burn;  this  material  is  ultimately  exhausted 
without  further  reaction.  However,  some  fraction  of 
these  gases  can  be  partially  burned,  resulting  in 
products  of  incomplete  combustion  and  especially 
CO. 

A  full  kinetic,  spatially  varying  model  [99]  can  ac¬ 
count  for  the  temperature  distribution  in  the  bound¬ 
ary  layer  and  crevices,  in  addition  to  the  main  bulk 
charge,  to  predict  unburned  hydrocarbon  and  CO 
emissions  from  this  engine.  As  a  result  of  this  varying 
temperature  distribution  in  the  fuel/air  mixture,  the 
overall  heat  release  is  spread  out  in  time  relative  to 
the  one-zone  homogeneous  treatment,  in  much  bet¬ 
ter  agreement  with  experimental  results,  as  shown 
in  Fig.  8  for  cases  in  which  the  bulk  region  and 
boundary  layer  were  defined  by  10  spatial  zones. 
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Fig.  9.  Ratios  of  experimental 
[101]  to  numerical  results  [98]  for 
the  main  combustion  parameters  of 
HCCI,  for  natural  gas  at  different 
boost  pressures,  using  the  10-zone 
spatial  model. 
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Fig.  10.  Computed  concentrations  of  basic  species  in 
parts  per  million  in  the  bulk  gases  in  HCCI  ignition,  with 
natural  gas  fuel  [98], 

This  more  refined  model  describes  many  experi¬ 
mental  features  of  HCCI  very  well  except  for  the 
predicted  CO  emissions,  as  shown  in  Fig.  9. 

The  kinetic  details  of  these  model  computations 
indicate  that  HCCI  ignition  is  controlled  by  H202 
decomposition.  Computed  variations  in  selected 
species  concentrations  are  summarized  in  Fig.  10,  in 
which  it  is  clear  that  H202  decomposition  occurs  at 
the  time  of  autoignition,  observed  experimentally 
[100]  to  occur  at  about  1°  after  TDC.  The  fuel/air 
mixture  follows  a  temperature  and  pressure  history 
very  similar  to  those  encountered  by  end  gases  in 
spark  ignition  engines  and  in  a  rapid  compression 
machine.  Modest  heat  release  occurs  at  lower  tem¬ 
peratures,  depending  on  the  cool  flame  reactivity  of 
the  specific  fuel  or  fuel  mixture,  resulting  in  variable 
amounts  of  early  temperature  increase.  Variation  in 
any  engine  parameter  that  gets  the  reactive  fuel/ air 
mixture  to  the  H202  decomposition  temperature  of 
about  1000  K  earlier  will  advance  ignition,  and  any¬ 
thing  that  delays  reaching  that  temperature  will  re¬ 
tard  ignition.  Increasing  the  intake  manifold  tem¬ 
perature  and  use  of  proignition  additives  such  as 
ozone  (see  Fig.  4)  or  ethyl-hexyl  nitrate  advance  ig¬ 
nition,  while  exhaust  gas  recirculation  and  addition 
of  knock  suppressants  retard  ignition. 

In  many  ways,  HCCI  combustion  is  very  simple. 
Low  and  intermediate  temperature  reaction  se¬ 
quences  process  the  fuel/ air  mixture  during  its  com¬ 
pression;  the  amount  of  low  temperature,  cool  flame 
heat  release  varies,  depending  on  the  composition 
of  the  fuel.  Ignition  occurs  at  the  temperature  where 
the  core  fuel/air  charge  reaches  the  H202  decom¬ 
position  temperature,  so  models  developed  for  en¬ 
gine  knock  and  rapid  compression  machines  are  di¬ 
rectly  applicable  to  HCCI  systems.  Low  NOx 
production  is  a  result  of  the  very  low  overall  fuel/ air 
equivalence  ratio;  the  low  equivalence  ratio  is  pos¬ 
sible  because  flame  propagation  is  not  required  in 
the  HCCI  engine  [103],  Combustion  is  spread  out 


in  time  because  the  extended  boundary  layer  reacts 
later  than  the  core  gases.  High  hydrocarbon  and  CO 
emissions  are  an  inevitable  result  of  the  low  bulk  gas 
temperatures  and  wall  boundary  layers,  and  it  is  very 
likely  that  postcombustion  exhaust  gas  treatments 
will  be  necessary  to  make  HCCI  combustion  viable. 
If  low  temperature  hydrocarbon  catalysts  were  to 
become  available,  HCCI  would  represent  an  impor¬ 
tant  competitor  to  conventional  spark  ignition  and 
diesel  engines. 


Summary 

The  role  of  chain  branching  in  determining  the 
onset  of  ignition  has  identified  two  distinct  chain¬ 
branching  mechanisms  affecting  hydrocarbon  igni¬ 
tion  in  practical  systems.  The  high  temperature 
(above  ~1200  K)  mechanism  is  important  in  shock 
tube  ignition  and  practical  systems  including  deto¬ 
nations  and  pulse  combustors.  These  are  primarily 
thermal  ignitions,  where  heat  release  increases  the 
temperature,  further  increasing  the  rate  of  ignition 
until  reactant  depletion  ends  the  ignition. 

The  second  chain-branching  regime  depends  on 
H202  decomposition  at  about  900-1000  K.  Experi¬ 
ments  in  rapid  compression  machines  and  flow  re¬ 
actors  and  corresponding  model  calculations  con¬ 
firm  the  importance  of  H202  decomposition  kinetics 
and  provide  insight  into  the  mechanisms  of  ignition. 
When  these  models  are  then  applied  to  practical  sys¬ 
tems,  especially  in  engine  environments,  it  is  clear 
that  the  same  kinetic  features  control  ignition  in  all 
of  these  systems.  Modifications  of  the  combustion 
environment  that  enable  a  system  to  reach  the  de¬ 
composition  temperature  at  earlier  times  advance  ig¬ 
nition,  while  modifications  that  delay  reaching  that 
temperature  retard  ignition. 

This  analysis  shows  that  autoignition  in  engines  is 
dominated  by  only  one  elementary  reaction,  II202 
decomposition.  Low  temperature  reactions  or  cool 
flames,  including  the  highly  chain  branched  alkyl- 
peroxy  radical  isomerization  kinetic  system,  simply 
advance  the  time  at  which  H2Oa  decomposition  is 
observed.  This  simplifies  the  analysis  of  combustion 
ignition  phenomena  and  focuses  kinetic  attention  on 
a  very  limited  family  of  reactions  that  control  au¬ 
toignition  phenomena  in  a  veiy  wide  range  of  prac¬ 
tical  systems. 
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INVITED  TOPICAL  REVIEW 


COMMENTS 


Christian  Eigenbrod,  / ARM,  University  of  Bremen, 
Germany.  You  were  talking  nicely,  and  pretty  imaginatively, 
about  staged  ignition  in  diesel  engines.  But,  as  with  putting 
a  droplet  of  water  into  wine,  in  diesel  engines  the  physical 
processes  of  vaporization  and  mixing  are  not  completely 
over  when  autoignition  happens;  the  process  must  be  taken 
into  account.  As  remaining  droplets  pin  the  temperature 
in  their  vicinity  close  to  their  boiling  temperature,  the  gra¬ 
dient  field  affects  the  chemical  process  leading  to  autoig¬ 
nition  in  a  much  more  complex  way  than  homogeneous 
gas-phase  shock  tube  conditions  can  unveil. 

Author’s  Reply.  The  limited  time  available  for  my  oral 
presentation  resulted  in  my  inability  to  describe  the  diesel 
ignition  problem  in  the  detail  it  deserves.  As  noted  by  Dr. 
Eigenbrod,  vaporization  and  mixing  make  this  a  very  com¬ 
plex  problem.  The  overall  problem  and  our  kinetic  mod¬ 
eling  of  it  are  described  in  considerable  detail  by  Flynn  et 
ah,  based  on  the  experimental  laser  diagnostics  work  of 
John  Dec  from  Sandia  National  Laboratories.  That  mod¬ 
eling  analysis  reinforces  the  point  being  made  in  the  lecture 
that  the  key  chain-branching  reaction  path  in  diesel  ignition 
is  the  same  intermediate  temperature  path  that  controls 
engine  knock  ignition,  homogeneous  charge  compression 
ignition  (HCCI),  and  rapid  compression  machine  ignition, 
the  decomposition  of  hydrogen  peroxide. 

A  full  model  of  diesel  ignition  would  certainly  have  to 
include  fluid  mechanics,  liquid  fuel  vaporization,  and  spa¬ 
tial  transport  of  heat  and  chemical  species.  This  type  of 
model  is  being  attempted  at  the  present  time  by  a  number 
of  researchers.  The  point  of  the  present  work  is  that  sim¬ 
plified  kinetic  models,  making  a  number  of  assumptions, 
are  able  to  establish  the  major  reaction  paths  that  lead  to 
ignition.  Quantitative,  predictive  models  will  require  more 
detail  in  their  description  of  these  physical  processes. 

One  major  point  of  the  paper  is  that  shock  tube  ignition 
is  the  idealized  form  of  high-temperature  ignition,  while 
diesel  ignition  is  an  applied  example  of  intermediate-tem¬ 
perature  ignition.  Therefore,  the  two  processes  are  con¬ 
trolled  by  different  kinetic  reaction  mechanisms. 

• 

Arvind  Atreya,  University  of  Michigan,  USA. 

1.  From  your  excellent  lecture,  it  appears  that  ignition  and 
extinction  are  two  sides  of  the  same  coin.  I  agree.  Could 
you  comment  on  this? 

2.  If  the  above  is  true,  there  seems  to  be  a  discrepancy  in 
temperatures  at  ignition  and  extinction.  Experiments  by 
Tsuji  and  coworkers  show  that  the  temperature  at  ex¬ 
tinction  for  hydrocarbon  flames  is  about  1450  K, 
whereas  your  ignition  temperatures  are  around  1 100  K. 
Is  it  possible  that  extinction  follows  the  high  T  branch? 


Author’s  Reply.  In  another  paper  at  this  same  symposium 
[1],  we  used  kinetic  modeling  to  reproduce  the  work  of 
Tsuji  et  al.  cited  by  Prof.  Atreya,  which  was  also  reported 
by  Law  and  Egolfopoulos.  The  works  of  Tsuji  et  al.  and 
Law  et  al.  deal  with  extinction  of  atmospheric  pressure  hy¬ 
drocarbon  flames  at  temperatures  of  about  1450  K.  We 
extended  these  calculations  to  engine  pressures  of  50-100 
bar  to  show  that  the  extinction  temperature  increases  to 
about  1900  K,  with  profound  implications  about  emissions 
of  NO,.  Our  modeling  work  attributes  this  extinction  to 
suppression  of  the  high-temperature  chain-branching  re¬ 
action  pathway,  as  Prof.  Atreya  indicates.  The  ignition  tem¬ 
peratures  of  1100  K  cited  in  his  question  describe  different 
problems  including  engine  knock  and  early  premixed  ig¬ 
nition  in  diesel  engines,  not  flame  propagation.  Therefore, 
I  disagree  with  his  proposition  that  ignition  and  extinction 
are  parts  of  the  same  problem. 
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Norbert  Peters,  RWTH  Aachen,  Germany.  In  the  inter¬ 
mediate  temperature  range  of  n-heptane  ignition,  radical 
runaway  at  the  second  ignition  is  triggered  by  the  depletion 
of  the  fuel  which  is  the  chain-breaking  agent  during  the 
period  after  the  first-stage  ignition.  Hydrogen  peroxide  de¬ 
composition  provides  the  radicals  needed  for  fuel  con¬ 
sumption,  but  its  temperature  dependence  is  not  a  deter¬ 
mining  factor  because  the  temperature  is  nearly  constant 
during  this  period.  If  it  was,  hydrogen  peroxide  would  de¬ 
crease  rather  then  increase  during  this  period. 


Russell  Lockett,  City  University,  UK  A  dynamical  sys¬ 
tems  analysis  of  methanol  autoignition  reveals  that  thermal 
runaway  begins  at  the  maximum  of  the  hydrogen  peroxide 
concentration  profile.  This  maximum  also  defines  a  bifur¬ 
cation  in  the  system.  Therefore,  autoignition  is  defined  by 
the  decomposition  of  hydrogen  peroxide.  I  would  expect 
this  to  be  true  of  the  thermal  runaway  in  two-stage  ignition 
as  well. 

Author’s  Reply.  In  answer  to  the  previous  two  com¬ 
ments,  it  is  tempting  to  try  to  simplify  these  kinetic  systems 
to  make  them  easy  to  understand.  A  linear  model  is  fre¬ 
quently  most  satisfying:  A  leads  to  B  leads  to  C.  In  these 
kinetic  problems,  several  processes  happen  together  to 
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cause  the  thermal  runaway  of  ignition,  and  it  is  not  clear 
that  any  of  them  can  be  identified  as  the  "cause.”  It  is  most 
useful  to  recall  the  different  processes  that  occur  and  the 
conditions  required  for  each  of  them.  Hydrogen  peroxide 
requires  a  temperature  of  1000—1100  K  to  decompose,  and 
its  decomposition  provides  lots  of  OH  radicals.  Hydrocar¬ 
bon  fuels  react  rapidly  with  radical  species  such  as  OH,  O, 
and  H.  Rates  of  reaction  with  these  radicals  are  much 
greater  with  hydrocarbons  than  with  species  such  as  CO 
and  H2,  which  produce  rapid  heat  release.  The  presence 
of  hydrocarbon  effectively  inhibits  these  heat  releasing  re¬ 
actions,  so  hydrocarbons  inhibit  ignition. 


John  Griffiths,  University  of  Leeds,  UK  I  recognize  the 
need  for  economy  of  discussion  in  a  wide-ranging  review, 
but  I  am  concerned  by  the  simplification  that  there  is  a 
"critical  decomposition  temperature”  for  the  H202  inter¬ 
mediate.  This  is  open  to  misinterpretation  because  the 
event  is  not  a  “criticality”  but  a  subtle  interaction  between 
the  amount  of  H202  formed  and  the  evolving  temperature 
of  the  reactants.  As  you  have  noted,  the  influence  of  tem¬ 
perature  is  dominant  because  the  activation  energy  for  de¬ 
composition  is  so  high.  The  chain-thermal  feedback  as  the 
decomposition  proceeds  leads  to  an  acceleration  that  may 
look  like  a  “critical  condition”  for  temperature  (and  even 
concentration,  maybe),  but  it  is  not  so.  There  are  relation¬ 
ships  also  to  the  duration  of  events  in  the  lead  up  to  the 
"proper  ignition”  through  the  “starting”  temperature  and 
the  rate  at  which  H202  is  generated.  The  higher  the  “ceil¬ 
ing  temperature”  for  the  R  +  02/R02  equilibrium,  the 


shorter  will  be  that  second  stage:  the  system  is  given  a 
reduced  handicap  (which  is  governed  by  the  R02  struc¬ 
ture)!  The  H202  formation  obviously  relies  on  H02,  but  I 
am  inclined  to  single  out  CHoO  as  an  important  molecular 
intermediate  in  this  stage  because  all  routes  for  its  oxida¬ 
tion  lead  to  H02.  The  relevance  of  CHoO  in  the  develop¬ 
ment  of  the  second  stage  was  discussed  by  Gibson  et  al. 
[1]. 
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Author’s  Reply.  The  analysis  included  in  the  present  pa¬ 
per  on  critical  temperatures  for  decomposition  of  H202 
was  taken  directly  from  the  work  of  Prof.  Griffith.  He  and 
I  have  communicated  many  times  on  these  issues.  This 
comment  is  another  example  of  the  fact  that  kinetic  sys¬ 
tems  are  extremely  complex,  even  when  all  the  rate  con¬ 
stants  are  known.  Professor  Griffith  and  his  collaborators 
have  studied  the  mathematical  issues  associated  with  ig¬ 
nition,  providing  important  insights  about  critical  phenom¬ 
ena.  The  identification  of  the  ceiling  temperature  as  an 
important  metric  for  this  type  of  problem  has  been  dis¬ 
cussed  at  length,  and  it  is  clearly  important.  The  idea  of 
the  importance  of  formaldehyde  as  a  major  molecular  in¬ 
termediate  is  interesting  but  not  yet  supported  by  any  ki¬ 
netic  evidence.  This  would  be  a  suitable  topic  for  future 
research. 


